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We discuss the dissipative diffusion-type term of the form m x V^ftm in the phenomenological 
Landau-Lifshitz équation of ferromagnetic precession, which describes enhanced Gilbert damping of 
finite-momentum spin waves. This term arises physically from itinerant-electron spin flows through 
a perturbed ferromagnetic configuration and can be understood to originate in the ferromagnetic 
spin pumping in the continuum limit. We develop a gênerai phenomenology as well as provide 
microscopic theory for the Stoner and s-d models of ferromagnetism, taking into account disorder 
and electron-electron scattering. The latter is manifested in our problem through the Coulomb drag 
between the spin bands. The spin diffusion is identified in terms of the transverse spin conductivity, 
in analogy with the Einstein relation in the kinetic theory. 
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I. INTRODUCTION 

The problem of spin diffusion through conducting fer- 
romagnetic médium attracted much attention over sev- 
eral décades. Semiclassical spin transport in the 
présence of a weak magnetic field H can be captured by 
the conventional diffusion équation (neglecting spin re- 
laxation) : 

ôtS = H X S + DV^S, (1) 

where D is the diffusion coefficient and H is the total 
effective field (omitting the gyromagnetic ratio), includ- 
ing the applied and exchange contributions. The first 
term on the right-hand side describes spin precession in 
the local field while the second term stands for the ordi- 
nary diffusion of spin density S. Equation ^ is, how- 
ever, not applicable to most realistic ferromagnets, whose 
spin interactions are characterized by a large exchange 
energy Axe- In particular, when Axe is comparable to 
the Fermi energy (which is the case in transition met- 
als), the spin precession in the exchange field cannot be 
treated in the diffusive transport framework. Further- 
more, the time- depe ndent exchange field induces spin- 
pumping current^^'*' inside the ferromagnet with spatially 
inhomogeneous magnetization dynamics, which can con- 
siderably modify the self-consistent magnetic équation of 
motion. Here, we wish to elucidate the central rôle of 
such self-consistent dissipative spin currents, which gov- 
ern the diffusion-like terms in the magnetic équation of 
motion in the limit of strong ferromagnetic exchange cor- 
relations. 

This paper is a foUow up to our previous work;^ pro- 
viding additional technical détails and offering a broader 
phenomenological base. Apart from assuming strong 
exchange corrélations limit, our phenomenological ap- 
proach and the main results of the paper should not be 
sensitive to the microscopic détails and do not rely on 
the spécifie model of the ferromagnetic material (such 



as the Stoner or an s-d model, for example). The main 
goals of this paper are as foUows: (i) to put the results of 
Réf. El into a broader phenomenological perspective, (ii) 
to explicitly show that two quite différent models — the 
spin-polarized itinérant électron liquid (treated in Réf. lïï]) 
and the s-d model — lead to the same phenomenology and 
can be treated in parallel, and (iii) to make direct contact 
with the spin-pumping theory.EEl 

To be spécifie, let us consider a continuons ferromag- 
netic médium, with the effective field and spin density 
initially pointing along the z axis. For weak excitations 
close to this state, we may try expanding the ensuing 
transverse spin-current density as^ ^ 

i, = -D'^xdS-D"dS, (2) 
which enters in the continuity équation: 

i—x,y,z 

In the limit of vanishing ferromagnetic corrélations, we 
recover Eq. ^ by setting D' ^ and D" ^ £) in Eq. g . 
Hereafter, we are focusing exclusively on the transverse 
spin dynamics and spin currents. The longitudinal spin 
flows are conventionally described in terms of the ordi- 
nary diffusion for spin-up and spin-down électrons with 
spin-dependent diffusion coefficients and spin-flip scatter- 
ing between the up and down spin bands.lïûlUnderstand- 
ing the transverse spin flows and dynamics requires more 
care, in part due to the inherently quantum-mechanical 
behavior in the case of a strong exchange field. When 
the magnetic excitation is driven by the self-consistent 
transverse field h = z x H x z, there should also be field- 
driven contributions to the transverse spin current Q, 
such as ji c>c dih.. 

The problem in fact simplifies in the limit of strong 
exchange corrélations. We will in the foUowing employ 
a mean-field view of ferromagnetism, where the collec- 
tive spin dynamics are driven by the exchange field, H = 



2 



— Axcïii(r,t) (setting h = 1 throughout), parametrized 
by the local and instantaneous spin-density orientation, 
m = S/ S, which has to be solved for self-consistently. 
Since we are only interested in the transverse spin dy- 
namics, we set the magnitude of the spin density S to be 
spatially and time independent. In the limit of large Axe, 
the spin currents can be parametrized by m(r, t). We can 
thus proceed phenomenologically and expand ji in spa- 
tial and time derivatives of m(r, t). For a static magnetic 
profile m(r), we have the familiar exchange spin flow 



j'j = —A m X â^m 



(4) 



(where A is the material-dependent exchange-stiffness 
constant), which is the only first-order form allowed by 
spin-rotational and time-reversal symmetries. To avoid 
unnecessary complications, we will assume isotropic fer- 
romagnet throughout this paper. Dynamics allow for 
dissipative spin-current contributions that break time- 
reversal symmetry: 



-7/ m X didtm ■ 



(5) 



Focusing on linear déviations of m from the equilibrium, 
m(°) = z, we omit ternis such as diina x ôtm. 

According to the time-reversal property, the spin- 
current density Q corresponds to the D' term in Eq. (|2| , 
while the spin-current density ([5| is analogous to the D" 
term, although the latter two are certainly not idcnti- 
cal. In fact, we wish to emphasize the striking différ- 
ence between the diffusive picture for the spin currents, 
Eq. ([2]), on one side and Eqs. Q and ^ on the other 
side, where we expand spin currents phenomenologically 
in terms of the time-dependent magnetic texture m(r, t). 
The latter approximation is spécifie to the limit of strong 
exchange corrélations, where the nondissipative spin cur- 
rent, Eq. Q, is determined by the instantaneous mag- 
netic profile, while the dissipative spin current, Eq. ([5|, 
can be interpreted as quasiparticle spin pumping by the 
collective magnetic dynamics,* rather than ordinary spin 
diffusion. It is also instructive to draw analogy between 
coefficients A and 77 in Eqs. Q, ([s]) and the shear modu- 
lus and shear viscosity, respectively, in elasticity theory. 

In the next section, we develop further the phenomeno- 
logical grounds for Eqs. Q and ([s]), before proceeding 
with microscopic calculations for the dissipative coeffi- 
cient 77 in Secs. III and IV In Sec. |V] we discuss a spin- 
pumping interprétation of dissipative spin current ([s]), 
before summarizing our work in Sec. |VI[ 



II. PHENOMENOLOGY 

A. Landau-Lifshitz theory 

The conventional starting point for studying ferro- 
magnetic precession is the nondissipative Landau-Lifshitz 
(LL) équation 



where we define the effective field H* as the functional 
derivative of the free energy: 



n* = d^F[m]/S. 



(7) 



In this Landau-Lifshitz phenomenology, which is applica- 
ble well below the Curie température, only the position- 
dependent direction of the magnetization is taken to 
be a dynamic variable, parametrizing the Free energy 
F[m(r)]. The angular-momentum density S = Sm is 
assumed to be related to the magnetization by a con- 
stant conversion factor, the effective gyromagnetic ratio. 
(Abusing terminology, we say spin density synonymously 
with angular-momentum density.) Since in the common 
transition-metal ferromagnets the gyromagnetic ratio is 
négative, we wrote Eq. ([t]) with an extra minus sign in 
comparison to the standard définition, where m is taken 
to be the direction of the magnetization rather than the 
spin density. The right-hand side of Eq. (|6| is the phe- 
nomenological reactive torque on the spatially-resolved 
magnetic precession, which generalizes the simple Larmor 
precession of Eq. Q. Note that the dissipation power 



P[m(r,t)] ^-S j (f'rW ■ dtm. 



(8) 



clearly vanishes according to Eq. ([ô]). We also eas- 
ily verify that the time reversai (under which t — > —t, 
m —m, and H* —H*) leaves Eq. (|6| unchanged, 
as it should in the absence of dissipation. The only dis- 
sipative term we can write in the quasistationary limit 
(i.e., up to the first order in dt), assuming spatially uni- 
forni and isotropic ferromagnet, is the so-called Gilbert 
daniping.'ï^ 



ami 



LLG 



H* X m — a m X dtm , 



(9) 



where a is a material-dependent dimensionless (Gilbert) 
constant. A typical expérimental value for a turns out 
to be often of the order of 10~^ in varions metallic ferro- 
magnets, which means that it takes roughly 2TT/a ^ 10 
precession cycles for an out-of-equilibrium magnetiza- 
tion to relax to a static equilibrium direction along H*. 
The Gilbert damping breaks time-reversal symmetry and 
causes a finite dissipation power: 



P[m(r,t)] = aS / (^m) 



ami 



LL 



H* X m, 



(6) 



(10) 



As a side comment, we note that an alternative, so-called 
Landau-Lifshitz damping term m x H* x m is mathemati- 
cally identical to the Gilbert damping mx^m in Eq. (|9|, 
up to an extra factor of (1 -|- a^) on the left-hand side of 
the équation. 

The effective field H* is in practice dominated by the 
applied magnetic field, magnetic crystal anisotropies, and 
magnetostatic (dipole-dipole) interactions. In the prés- 
ence of spatial inhomogeneities, there is also exchange 
contribution to the free energy, which to the leading 



3 



(quadratic) order in magnetic inhomogeneities can be 
written as^-'^' 

P^c = ^ I d\ [{d^mf + {dyxnf + (9,m)2] . (11) 



The corresponding effective field is 

H,c = -(v4/5)V2m, 
and the associated terni in LL Eq. ^ is 

dtm\^^ = (A/5) m X V^m. 
This équation can also be formally written as 



(12) 
(13) 



i—x,y,z 

— ~A m X diTcn , 



(14) 



which simply recovers our equilibrium spin current Q. 
We emphasize that tfiis spin current does not dépend on 
magnetic dynamics. 

To summarize thèse prehminary considérations, the 
phenomenological LL équation describes collective mag- 
netic precession driven by local effective fields as well as 
equilibrium spin currents. At this point, there is, how- 
ever, a conspicuous asymmetry in the treatment of the 
dissipative correction to the LL équation, i.e., Gilbert 
damping ([Ô]), which dépends only on the local magnetic 
dynamics and thus does not involve spin currents. To 
overcome this "discrepancy," we expand the dissipative 
terms to second order in spatial derivatives, generalizing 
Gilbert term to 



-a m X 9tm + (îz/S') m X V^ôtm , (15) 



where ?7 is a new phenomenological parameter, charac- 
terizing spin-wave damping. Assuming spatial-inversion 
symmetry (under which dt — > —di and m — > m), pre- 
vents us from writing any phenomenological terms linear 
in spatial derivatives. Recall also that we are always as- 
suming small perturbations with respect to a uniform 
equilibrium magnetization, so that ail spatial and time 
derivatives must hit a single m (for example, a dissipa- 
tive term of the form J^ii'^i^ ' (™ ^ dtm.)]dim is dis- 
regarded since it is higher order in small déviations of 
m). Additional quadratic terms would be allowed phe- 
nomenologically if, e.g., we developed our linearized the- 
ory with respect to an equilibrium magnetic texture, such 
as a domain wall or magnetic spiral. Some of such terms 
were discussed in Réf. 13, which is beyond our présent 
scope. Finally, we note that we wrote Eq. (15) with no 
direct coupling to the effective field H* . We justify this 
by assuming that the ferromagnetic corrélations are char- 
acterized by a very large energy scale Axe, so that mi- 
croscopic processes responsible for dissipation are driven 
by the collective variable m, rather than directly by H*. 
In transition- métal ferromagnets, the internai exchange 
energy is of the order of eV, while the effective field H* 



corresponds to microwave frequencies (i.e., at least three 
orders of magnitude smaller than the exchange energy) . 
This means that when we excite magnetic dynamics by 
an external field, the microscopic degrees of freedom re- 
spond not to the small driving field but rather the much 
larger self-consistent exchange field parametrized by the 
time-dependent m. For the same reason, the spin current 



in Eq. (14) dépends only on the magnetic profile m(r). 



irrespective of how it is created by applied fields. 

The total dissipation power corresponding to Eq. (15) 
now becomes 



P[m(r, t)] = J (fr [aS [dtraf + 77 {^^^traf 



(16) 



Similarly to Eq. (14), we can also write the 77 term in 



Eq. ( 15 ) in the form of the divergence of the spin-current 



density 



-77 di{m X dtm) , 



(17) 



reproducing Eq. ([s]). We thus identified two contribu- 
tions to the spin-current density: usual exchange spin 
current ( 14 1 and dissipative spin current ( 17 1 , which 



we will late r interpret as the dynamically-driven spin 
pumping.l^ Spin current (17) can thus damp down 



spin-wave excitations even in the absence of any spin- 
relaxation scattering.t23l The latter is, however, believed 
to be the culprit for a finite Gilbert damping a,^^ which 
relaxes uniform magnetic precession by transferring its 
angular momentum to the atomic lattice. 



In the présence of dissipative currents ( 17 ) , the relative 
linewidth of the spin-wave resonanc^ï^ is proportional to 
a + {ri/S)q^, for the wave vector q. In the absence of 
the Gilbert damping a, thus, the spectral width of the 
spin-wave excitation would vanish in the long-wavelength 
limit.'ij 



B. Mermin ansatz for spin current 

We now wish to establish a microscopic procédure for 
eval uati ng the dissipative component of the spin current, 
Eq. (17). Ref . [ni adapted Mermin ansati^^ for this pur- 
pose, which we will reproduce below. 



Microscopically, 
the spin-current density ji is carried by conducting élec- 
trons responding to the mean-field exchange interaction 



-Axcm(r,t) • à/2 



(18) 



in the self-consistent single-electron Hamiltonian (which 
could stem, e.g., either from the coupling to the localized 
d électrons in the s — d model or the itinérant électron 
Stoner/LDA exchange) . <t is the vector of Pauli matrices, 
which defines the électron spin operator. 



Let us for the moment view exchange interaction ( 18 1 



as an external parametric driving field, not concerning 
with a self-consistent détermination of m(r,i). In par- 
ticular, we may allow for an instantaneous déviation of 
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the électron spin density s from the exchange-field direc- 
tion m. This will allow us for a trick to find the ensuing 
spin flows, which is what we are after. The spin-density 
continuity équation corresponding to Hamiltonian ( 18 1 is 



dts = Axe z X (s m - s) - diji 



(19) 



The equilibrium orientation of m is taken to be along the 
z axis and we assume small-angle excitations, which do 
not modulate the magnitude of the spin density, s = |s|. 
s here is the spin density of the conducting électrons, 
which in, e.g., the s — d model has to be distinguished 
from the total spin density S that enters Eq. 

We next use the Mermin ansatz to relate the spin- 
current density ji to the spin density s: 



j., = a±A^cdi{m - s/s) 



(20) 



where a± is the transverse spin conductivity, to be evalu- 
ated later by the Kubo formula. Eq. ( 20 ) is analogous to 



Ohm's law for electric current density, with the expres- 
sion on the right-hand side reminiscent of the gradient 
of the electrochemical potential. The physical reasoning 
behind ansatz (20l is simple: there should be no dis- 



sipative spin currents in the static configuration, which 
corresponds to s(r) = sm(r). The advantage in writing 
the spin current in this form is that a± will now have 
to be evaluated in the limit of (q, w) 0. Combining 
Eqs. (19l and (20l will then give us the spin current to 



the linear order in q and ui: exactly what we need to re- 
late it to Eq. ([T7| and read out 7]. In fact, it is sufïicient 
to find Axe (m — s/s) « — z x dtm from Eq. ( [T9| , which 
is valid to the linear order in lu and zeroth order in q, 
before putting it into Eq. pÔ] ) to finally find 

(21) 



ji = -CT_L di{z. X dtm) . 



Comparing this with Eq. ( 17l, we immediately identify rj 



with the transverse spin conductivity: 
7] = cr_L . 



(22) 



Equation ( 22 ) can be interpreted as an analog of the Ein- 



stein relation for transverse spin diffusion in strong fer- 
romagnets. 

C. Transverse spin conductivity 

As is the case with the charge conductivity, it is con- 
venient to evaluate the transverse spin conductivity in 
the velocity gauge. Namely, we eliminate the spin "po- 
tential," corresponding to small magnetization déviations 
ôm = m — z in Eq. (18), by the SU(2) gauge transfor- 
mation 



V'(r,t) 



JA„c Jl^dt' Sm(r,t')-à/2 J r 



V;'(r,i) 



(23) 



at the expense of introducing the SU(2) vector potential 



A,, = -A, 



dt'9im(r,t') • <t/2, 



(24) 



which enters the kinetic part of the single-particle Hamil- 
tonian as 



(25) 



where Pi = —idi and m* is the electron's effective mass 
(assuming exchange-split parabolic bands). It is easy to 
verify that the effe ctive field driving the spin current in 
velocity gauge (25), Ei — —dtAi, is the same as the fic- 
titious field Ei = —diV in original length gauge (18). 



One caveat is in order: Eqs. (23l-(25l are only valid for 



an Abelian exchange potential, which would be the case 
if only one vector component of (5m(r, t) was modulated 
(e.g., Srrix or 5my) in space and time. Such scénario is 
sufhcient for our purpose, in order to establish the trans- 
verse spin conductivity entering Eq. ( 20 ) . 

Fourier transforming the electric field Ei in time, 



the usual relationship is obtained: Ei{u}) 



iuiAi{uj). We now proceed to construct the semiclassi- 
cal transport équation for the spin current driven by a 
spatially homogeneous fictitious field = Tr[ijiô'] = 
Axcôiin, to deduce the long-wavelength conductivity de- 
fined by Ohm's laW^ 



(26) 



The semiclassical spin-current response, in the présence 
of the exchange splitting Axe, with disorder and electron- 
electron scattering is given bjŒ^ 



dtii + AxcZ X = 



4m* 



1 

-dis 



(27) 



where n is the total equilibrium (conducting) électron 
density. The second term on the right-hand side of 
Eq. ( 27 1 describes spin-current relaxation, due to dis- 
order and electron-electron scattering. Note that even 
in Galilean-invariant Systems, spin-independent Coulomb 
interaction between électrons causes relaxation of a ho- 
mogeneous spin current, in contrast to the ordinary cur- 
rent. Solving Eq. (27) at low frequencies, we recover 



Eq. ( 26 ) for the current component along , with^ 



4m* 1 + (TiAxe)2 



(28) 



where the total transverse spin scattering rate is defined 

by 



1 



1 

r-dis 



1 



(29) 



In particular, in the limit of weak spin polarization and 
no electron-electron interactions, t_l should reduce to the 
ordinary momentum scattering time t, and ctj^ to the 
quarter of the Drude conductivity nr/m* . 
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III. MICROSCOPIC CALCULATION 



B. Spin- force autocorrelator 



A. Spin-current autocorrelator 

In order to substantiate the preceding phenomenology, 
we need to establish the microscopic expressions for the 
involved scattering times, r^'^ and t'^. In the velocity 
gauge discussed in the previous section, the transverse 
spin conductivity is given, according to the Kubo for- 
mula, by the spin-current autocorrélation functionP 

where the summation is over ail électrons in volume V 
and 

/•OO 

{{Â-B))^ = -i dte^(-+«+)*([i(t),S(0)]) (31) 
Jo 

represents the Fourier-transformed retarded (Kubo) 
linear-response function for the expectation value of the 
observable A under the action of a classical field that 
couples linearly to the observable B. 3to in Eq. (30l is 
inserted out of convenience, since the linear in lo response 
function is guaranteed to be imaginary. (The zeroth- 
order in w correlator includes also the omitted "diamag- 
netic pièce" of the spin current in the velocity gauge.) 

Assuming isotropic disorder (and for the moment no 
electron-electron interactions), the ladder vertex correc- 
tions to the conductivity vanish and we only need to eval- 
uate the bubble diagram defined by the (single-particle) 
spin-dependent Green's functions 



In the présence of electron-electron interactions, it 
is convenient to express the s^nn- current autocorrelator 
( |30[ ) in terms of the spin-/orce autocorrelator. To this 
end, we use the équation of motion for the operators 



defining Kubo formula ( 30 1 to find 



^((Si ^xlFxl] El Ô'xlFxl))u 



— ^1 o ■ T -.^ lim ■ 

(35) 

where Fxi = Pxi — —'i-[Pxh H] is the force operator along 
the X axis for the Ith électron. Evaluated with respect to 
a uniform magnetization, m = z, the force operator Fxi 
consists of two pièces: the disorder force and the electron- 
electron interaction force. Evaluating correlator ( 35 ) in 



the clean limit to second order in Coulomb interactions, 
one finds for the transverse spin scattering ratel^ 



1 



T{p)m*alri{kBTy^ 



(36) 



where is the Bohr radius, T température, kg Boltz- 
mann constant, the dimensionless Wigner-Seitz radius, 
and T(p) is a dimensionless function of the spin polariza- 
tion p = (n| — n^)/n {rig being spin-s électron density), 
which was discussed in Refs. I5|9I Notice that scattering 



rate ( 36 1 has the Landau quasiparticle scaling with tem- 
pérature. The finite-frequency modification of scatter- 
ing rate (|36| is, furthermore, accomplished by replacing 
(27rfcBT)^ {2'KkBTY + cj2. 



C. Spin-density autocorrelator 



G^'^(P,^) = 



1 



(32) 



where a =t, i (= ±) is the spin index along the z axis, 
Pfj = /z-|-crAxc/2 is the spin-a électron Fermi energy, is 
the chemical potential, and t„ is the spin-dependent dis- 
order scattering time. In the Born approximation for 
dilute white-noise disorder, the scattering rate is pro- 
portional to the électron density of states, and we can 
Write To- — Tvjvc, where t parametrizes the strength of 
the scattering potential, v„ is the spin-cr band density of 
states, and i' = {im + i'i)/2. A straightforward calcula- 
tion then leads tcP 



crj_ = 



4m* rf "A2^ ' 



in the strong exchange coupling limit, where 



1 

-dis 



(33) 



(34) 



3 nr(i/| ^ + ly^^) 
identifies the disorder contribution to effective transverse 



spin scattering rate ( 29 1 



It is also possible to calculate the transverse spin dif- 
fusion directly, as a linear spin-density response to the 
transverse magnetic field. We will carry that out in 
Sec. |IV| for two popular mean-field models of ferromag- 
netism in metals: the Stoner and the s — d models. In ad- 
dition to offering an alternative approach to the problem, 
this dérivation provides a justification for the preceding 
heuristic utilization of the Mermin ansatz. 

Starting with the mean-field Hamiltonian for itinérant 
électrons 



H ^ 



P 

2m* 



i7(r)-M-A,e àj2, 



(37) 



and directly solving for the self-consistent spin-density 
response to a small driving magnetic field, we will dérive 
in the next section the foUowing gênerai relation: 



Axe r 3TOX+-(q,t^) 
rj — — ^ lim 



(38) 



valid at long wavelengths, ç — > 0. The axially-symmetric 
(Kubo) spin-response function is defined by 

X+-(q,^) = -^((s+(r,i);s_(r',0)))q,<., (39) 
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where s± = ± iSy is the transverse spin density of where v = {ly^ + i'i)/2 is the spin-averaged density of 
itinérant électrons. The disorder potential U (r) entering states at the Fernii level and r is the characteristic scat- 
Eq. (37l is, as before, taken to obey the Gaussian white- tering time. 



noise corrélations: 



(C/(r)f/(r')) 



1 



27r i/r 



(40) 



Writing the spin density s(r) = Tr [ô-/5(r)] /2 in ternis of the électron density matrix pQ^(r) = 'î'J^{r)'i> a{r) in spin 
space, we proceed to evaluate in the standard imaginary-tinie formalism. At température T, we have: 



X+-(q, i^n) 



T 
2V 



X! ^i(p + p' + ï'^m + î^n)^T(p'i p; *^™) 



pp ,m 



where 



g,{p^p';^co^) = -^ J dWr'^'^'^dre-'P"-+'P''^'+''-''-^(vI/,(r,r)vI/t(r',0)) 



(41) 



(42) 



is the finite-temperature single-particle Matsubara Green's function. il„ = 2mTT is the bosonic and ujm = (2m + l)7rT 
ferniionic Matsubara frequencies, where n and m are integer indices. 
The disorder-averaged Green's function is given by 



lujjn - £pa + isign(w„)/2r<j 



(43) 



where Spa- — /2m* — /i — crAxc/2. The analytic continuation of the Matsubara Green's functions into the retarded 
(advanced) Green's functions is accomplished by replacing icu^ — > a; ± iO'^ and sign(ci;„j) — > ±. According to our 
convention (40l, t„ = Tvjv^. Taking into account the vertex ladder corrections (as shown in Fig. [T|), we obtain for 
the disorder-averaged response function: 



Ep^i(P + q;«^m +îî^n)ST(P;îW„) 



(44) 



where Ç — 1/2ttvt and by the Green's functions with a single wave-vector argument here we understand disorder- 
averaged propagators (43l. Inserting Eq. (43 1 into Eq. (44l and performing an analytic continuation onto the real 



frequencies, it is straightforward to calculate {(\,uj). Setting the température to zéro and taking the w limit, 

we find: 



STOX+_(q,w) 



A A 



47r (l_^x^A)(l_^^AA) 



(45) 



where {q) — J dpG'^ {p)Gj {p — q) and J dp = J d'^p/{2TT)^ in three dimensions. AU énergies entering thèse 
Green's functions are set at the Fermi level. To the lowest order in l/rA, we now obtain: 



3m x+-(q,t^) 



UJ 

8tt 



J dpA^{p)A^{p~q)+i^Qm J dpGfip-q)A^{p) J dpGf{p-q) gîeGf(p) 



where Ag- — —2^mG^ is the spectral function. 



r 



(46) 



The second term in Eq. ( |46| is the vertex ladder correc- 
tion, which is necessary for Eq. ( 46 1 to give a meaningful 



resuit. In particular, the vertex correction cancels the 
spurious g = contribution of the first term, which would 
give a 1/rAxc. Finally, in the limit of ç ^ ttAxc/wf, 



we arrive at: 

Q'"T,x+-(q,cj) 



3m*T(i.-i + i.;i)A4^'^'^ 



Using Eq. (38l, this finally gives: 



V 



(47) 



(48) 
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which agrées with Eqs. ( 22 1 . ( 28 1 , and ( 34 ) in the relevant 
here limit of rj^^ ^ Axe- 



IV. MEAN-FIELD FERROMAGNETISM 

A. Time-dependent LDA 

In a spin-density-functional thcory (s-DFT)l^li^ the 
many-body problcm of itinérant ferromagnetism is re- 
duced to the single-electron Hamiltonian 



Hit) 



2m' 



+ f/(r) - ^l 



[Axe m(r, + cjQZ + h(r, t)] • <t/2 . (49) 



ujQ <^ Axe is the ferromagnetic Larmor precession fre- 
quency in the présence of a uniform magnetic field ap- 
phed along the z axis. Axe m(r, t) is the self-consistent 
exchange field, such that Hamiltonian (49) produces the 



correct spin-density response. Since we are ultimately in- 
terested in the équation of motion for the collective ferro- 
magnetic dynamics, the spin-density response is ail that 
is needed. In the local-density approximation (LDA) of 
the s-DFT, the exchange field foUows the local and in- 
stantaneous magnetization direction m{r,t). h.{r,t) is 
the external rf driving field, which we will treat pertur- 
batively. 

The time-dependent portion of the Hamiltonian is thus 
given by 

È'{t) = - [ (fr [Axe Sïn{r, t) + h(r, t)] ■ à/2 . (50) 



Since ôm = Ss/S (where 6 dénotes small déviations from 
equilibrium) , we have for the transverse spin component 



s+(q,w) = x+-(q,t^) 



h+{q,uj) + -|^s+(q,w) 



(51) 



The self-consistent response function to the rf field, 
Xa — s+lh+, is thus 



-1 / \ ~-l / \ '-^xe 

x+-(q,'^) = x+-(q,t^) - -g- 



(52) 



In the LDA approximation, the problem thus trivially 
reduces to calculating the spin-spin response function for 
a noninteracting Hamiltonian with a fixed exchange field. 
Let us in gênerai write 



x+-(q,w) 



s 



[wr(q, ^) + Axe ia(q, lo)uj 



(53) 



in terms of functions and a that are to be determined. 
The self-consistent response function then becomes: 



x+-(q,w) 



s 



[wr(q, uj) — Lo\ — iQ!(q, lû)uj 



(54) 



At q = 0, obviously a;r(w) = wo and a(w) = 0. This fol- 
lows in gênerai from the spin conservation in the présence 
of Coulomb interactions and arbitrary spin-independent 
potential U{v). In this paper, we are most interested in 
the q-dependent damping function Q;(q, lo), which can be 

identified by a microscopic évaluation of (qi'^)- 

the limit of strong exchange corrélations. Axe S> w^, we 



immediately obtain from Eq. ( 53 1 



a(q,w -> 0) 



^lim^ 
S 



(55) 



In inversion-symmetric Systems, the leading in q spin- 
wave contribution to Gilbert damping is a(q, w ) = 
{rj/S)q^, so that self-consistent response function (54l 
corresponds to the dissipative term 

a^mldi,, = {i^/S)m X V^a*m (56) 

in Landau-Lifshitz Eq. ([6]) of motion for the magnetic 
spin direction m(r,É). This is the désirable resuit and, 
according to Eq. ( 55 1 , the microscopic expression for 77 



gives Eq. ( 38 ) of the previous section. In the next section. 



we will demonstrate that Eq. ( 55 ) is generic to mean-field 



treatment of conducting ferromagnets. 



B. s d model in RPA 

It is also instructive to pursue a more basic descrip- 
tion starting with a ferromagnetic lattice of localized d 
électrons exchange-coupled to itinérant s électrons. The 
corresponding Hamiltonian is 

È{t) ^Èo~J2 ["^•^^ • ^) + ■ ^)] ' (^^) 



where «Si are local d spins and Hq consists of the de- 
coupled Hamiltonian for itinérant électrons, de Zeeman 
Hamiltonian of the d électrons, as well as the d — d ex- 
change and possible dipolar interactions, h is the applied 
rf field, which we take for simplicity to couple to the lo- 
calized spin only. As long as the average exchange field 
experienced by the s électrons is sufficiently strong and 
the magnetization is dominated by the d électrons, we 
can disregard their direct rf coupling for our purpose. 
If also the Fermi wavelength is long in comparison to 
the d lattice spacing, we will treat the electronic band 
structure in the effective-mass approximation, and also 
coarse grain the local spins: J^i^-i ~^ I d^rS{r) and 
Si {V/N)S{r), where N/V is the density of d sites. 

Let us compute the spin-density response function for 
the d lattice: 

X+-(r, r'; t) = -^((5+(r, i); 5_(r', 0))) . (58) 

For this purpose, it is convenient to define bosonic 
magnon operators: 



(59) 
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q, i^n q, i^n 



î îlx'r, 




P + q ^ p + q 



FIG. 1: In this RFA summation for the magnon propagator, Eq. (64l, the magnon self-energy is provided by the s-electron 
spin-response bubble, Eq. |41l. Straight double lines dénote disorder-averaged s-electron Green's functions and dashed lines 
describe vertex ladder corrections. Each bubble contains ladder corrections to the vertex function. Wavy lines are the d-electron 
spin propagators. 



which obey the canonical commutation relations, 
[opiOp,] — (5pp', close to the fuUy-magnetized ground 
state. To be spécifie, let us take the Heisenberg model for 
exchange coupling, so that in the ground state, Si^ = T), 
the d-orbital spin [assuming the applied de magnetic field 
to point along the — z direction, as in Eq. (49l]. The d- 



orbital Hamiltonian for magnon excitations close to the 
ground state can thus be written as: 



(60) 



In terms of the magnon operators, we, furthermore, 
rewrite the s — d exchange interaction as 



J 



VN ^ 



y V 2 ^ [«P«+(P) + apS-(-p) 
p 

- JY,S^^S,{Vi) , (61) 



where s±(p) = / d'^r e^*P ""s±(r) is the Fourier- 
transformed transverse s-electron spin density. Approx- 
imat ing Siz ~ î?, we can combine the second term in 



Eq. (61| with iïo for the s électrons, resulting in the 
longitudinal mean-field exchange field Axe = JS, where 
S = VN/V is the averaged d-electron spin density. The 
spin-flop term of Eq. (61) can be expanded in terms of 



the electronic field operators as: 
J fVN 



H' 



V 



5][at.*|(p')*i(p' + p) 
pp' 



apvl/ (p')fT(p'-p) 



(62) 



In the imaginary-time (Matsubara) formalism, re- 
sponse function ( |58[ ) can be written as: 

(63) 



X+-(q,T) = -5" J'(q,T) , 
in terms of the magnon Green's function 
^(q,r) = -r(aq(T)at(0)) , 



(64) 



using définition (59l. Here, T is the time-ordering sym- 
bol. For noninteracting Hamiltonian (60), 

.Fo(q,ir!„) = — ^ , (65) 



where :P(q, if)„) = dre*""^J«^(q, t) and f]„ = 2n7rr 
as previously, is the bosonic Matsubara frequency. In or- 
der to calculate T in the présence of s — d exchange ( 62 1 



we will sum up the bubble diagramsfwhich constitute an 
RFA approximation) shown in Fig. [î] It is easy to recog- 
nize that each bubble (which is the magnon self-energy 
in this approximation) is just the itinérant électron spin- 
density response function that was already calculated in 
Sec. lIITCl 

Summing up the diagrams in Fig. [ï] we thus obtain 



T ^(q,ifi„) = \q, ifi„) - S(q,îî7„) 



(66) 



where the self-energy (the s-electron spin-response func- 
tion) 



S(q,zr!„) = -j25x+-(q,ïf^n) 



(67) 



foUows from Eq. (44l. Combining Eqs. (|63[), (66), and 



( |67| , we finally obtain (after analytic continuation onto 
real frequencies) : 



-_(q,-) = — 



-(q,t^) . 



(68) 



Equation (68) is actually quite trivial: it can be also ob- 



tained by treating the d-orbital magnetization dynamics 
and the associated s—d torque in a mean-field approxima- 
tion analogous to the preceding discussion of the Stoner 
model. Inverting Eq. (68), we can write it in form (54), 
after identifying 



w,.(q,w) = £q - S'J^5Rex-H-(q,w) 



and 



a(q, Lo) = SJ 



2 3?mx+-(q,w) 



iû 



(69) 



(70) 



which is, in fact, exactly the same as Eq. (55) in the 
low-frequency limit, using Axe = JS. We should em- 
phasize that, although is defined in this section for 



a différent physical System than xh in Sec. |IV A (and 

thus, not surprisingly, is found to be somewhat différent), 

the itinerant-electron response is the same through- 

out the paper. The reason why the magnetic damping 
a(q, oS) is identified in terms of the same quantity x^ 
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the two différent models of ferromagnetism can be traced 
to our phenomenological identification of this damping in 
ternis of the conducting-electron transverse conductivity, 
Eq. ( 22 1 . The latter is governed by the mean-field struc- 



b=L+d 



ture of the exchange field, irrespective of the microscopic 
origin of the ferromagnetic order. 

Let us also note in the passing that, unhke the ideahzed 
Stoner model considered in the previous section, the s — d 
magnetic damping may have a finite q = value even in 
the absence of any additional spin-dependent terms in the 
Hamiltonian. When the gyromagnetic ratios of the two 
électron species difïer (ultimately stemming from some 
form of spin-orbit interaction), the total spin no longer 
precesses undamped in the uniform field, and the uniform 
transverse spin component can decohere in the présence 
of ordinary scalar disorder. Since Eq. ( 70 ) corresponds 



to the magnetic field coupled to the d électrons only, we 
implicitly set the s électron g factor to zéro. 



V. SPIN-PUMPING INTERPRETATION 

It is illuminating to interpret the key resuit of this 
paper for the transverse spin diffusion of form (17) in 
terms of the spin pumping associated with a nonuni- 
form magnetic dynamics in ferromagnetic bulk.^ The 



ferromagnetic spin pumping was originally proposed in 
the context of magnetic multilayers with sharp normal- 
métal | ferromagnetic interfaces. This paper shows that 
analogous processes also take place in the continuons fer- 
romagnetic médium. 

To illustrate the direct connection between the trans- 
verse spin diffusion and the spin pumping, we consider 
a periodic stack of alternating F and N layers forming a 
two-component superlattice in the x direction.!^ We treat 
the model depicted in Fig. [2] in which an F|N bilayer 
forms the unit cell with thickness b = L + d, where the 
normal-metal spacer of width L séparâtes the magnetic 
films of thickness 



(71) 



The latter approximation allows us to neglect the trans- 
verse spin-current cohérence between two interfaces of 
the same magnetic layer.^ Translational invariance is as- 
sumed for simplicity in the latéral directions. We con- 
sider here collective spin-wave excitations, taking both 
the static and dynamic exchange couplings into account.l^ 
The static (RKKY-like) exchange interaction between 
neighboring ferromagnetic layers is mediated by the dis- 
sipationless spin currents flowing through the normal- 
metal spacer.ES We will parametrize the strength of this 
coupling by the corresponding precession frequency lu-^c 
of a single ferromagnetic film that is exchange-coupled to 
a pinned film. In the présence of the magnetic dynam- 
ics, additional dissipative spin currents set in. Their ori- 
gin lies in the spin pumping by the individual magnetic 
layers into the adjacent normal spacers, which at low 
frequencies is given bjEl ipump _ (?i/47r)gjj|pm x dtm. 



t 



t 



FIG. 2: A schematic view of the superlattice considered in the 
text: an F|N bilayer is repeated along the x axis, with either 
ferromagnetic or antiferromagnetic alignment of the consécu- 
tive magnetic layers. The System is translationally invariant 
along the two remaining axes. 



gjjjp is the dimensionless spin-mixing conductance per 

unit area of the F|N interface (which is assumed to be 
real-valued, for simplicity). This interfacial spin pump- 
ing induces nonlocal spin transfer in magnetoelectronic 
circuits, which can in gênerai be treated as a source term 
entering spin transport équations in normal and mag- 
netic layers. In a coUinear superlattice of Fig. [2] the 
problem simplifies considerably, because the spin-current 
vector IP""ip ex m X dtm is transverse with respect to 
the magnetic alignment (in both ferromagnetic and anti- 
ferromagnetic cases), within the linear-response régime. 
This means that the spin current pumped by one fer- 
romagnetic layer is either scattered back by the normal 
spacer and reabsorbed, or transmitted and absorbed by a 
neighboring layer, with no possibility to reach more dis- 
tant neighbors, subject to condition (71 ). Spin relaxation 



in normal spacers would only cause an overall increase 
in the effective Gilbert damping parameter of a uniform 
magnetic precession, and will thus be omitted, since our 
primary interest here is nonlocal damping effects. The 
problem of dynamic exchange between two adjacent fer- 
romagnetic layers thus effectively reduces to the analo- 
gous effect in magnetic bilayers, which was studied in dé- 
tail in Ref.|7] In particular, the net spin pumping through 
a given normal spacer is oc mi xôtmi — m2 xc}tm2, which 
reflects the dynamic spin injection in the opposite direc- 
tions by the adjacent magnetic layers mi and m2. Notice 
that the total pumping vanishes in the case of a perfectly 
synchronous precession, mi{t) = m.2{t). 

Let us now put the static and dynamic exchange inter- 
actions into the équation of motion for small-angle spin 
dynamics of a multilayer with respect to an all-parallel 
configuration, Ui{t) = mi{t) — z. For long-wavelength ex- 
citations, it may be approximated as a continuons func- 
tion u(a;, t) of the coordinate x normal to the interfaces. 
For the uniaxial effective field H* = —loq'l, the spin-wave 
dynamics obey the differential équation 

ôtu = [wqU - w^c&^ô^u -h aStU - ol\?dldtv\ x z , (72) 

where we made the foUowing définitions: Wxc = Jyic/ Sd 
and a' — 3fîe^p|j^|p/47rS'd. Here, Jxc is the exchange 
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coupling between two consécutive magnetic layers and 



(73) 



is the effective pumping résistance of the spacer. The first 
terni on the right-hand side of Eq. ( 73 1 parametrizes the 



pumping strength of the individual interfaces, as was dis- 
cussed above, and the second term is the ordinary ohmic 
résistance of the normal spacer (neglecting any spin relax- 
ation) , which backscatters the pumped currents and thus 
suppresses the dynamic exchange. The second spatial 
derivatives in Eq. ( 72 1 reflect simply the différence of the 



static and dynamic exchange spin currents through two 
consécutive normal spacers (which themselves require a 
finite misalignment of the adjacent magnetic layers) in 
the continuum limit. The static Heisenberg coupling can 
be interpreted as the superlattice équivalent of the bulk 
exchange-stiffness parameter A of Eq. ( 11 ), which for the 
superlattice becomes A — J^Jd- Both and a' are sen- 
sitive to the normal-interlayer thickness L, vanishing in 
the limit L ^ oo. It foUows from Eq. (72) that the 
small-momentum, q ^ 5^^, spin- wave excitations of the 
superlattice, propagating perpendicular to the interfaces, 
u ex exp[î((7a; — wt)], obey the dispersion relation 



l + i[a+ {bqYa'] 



(74) 



When ç ^ 0, oj{q) reduces to the Larmor frequency 
of the individual magnetic layers because the static and 
dynamic exchange couplings vanish when the consécutive 
magnetic layers move coherently in phase. Equation (74) 
holds up to momenta comparable to 6^^, when bq has to 
be replaced by 2sm{bq/2). 

The matters are quite différent for an 
antiferromagnetically-aligned superlattice, which is 
the ground state when, for example, Jxc < and 
H* — 0. In this case, we have a more complex 
dispersion: 



±^{bqY{l -f «2) _ 4^2 + j[2a + {bqfa'] 

~ ""''^^ l + a^+Aaa' + {bqYa'^ ' 

(75) 

where plus and minus signs refer, respectively, to the 
modes with antisymmetric and symmetric dynamics in 
the adjacent layers for overdamped motion, and to the 
right- and left-propagating modes when the real part 
of u){q) is significant. Note that now < 0, so that 
3m w > 0, as required for a stable configuration. In the 
absence of bulk magnetization damping, a = 0, Eq. (75) 
reduces to 



uj{q) 



±{bq)uj^c 
l±i{bq)a' ' 



(76) 



with linear dispersion and damping at small q. Equa- 
tions ( 75 1 and ( 76 ) can also be generalized to large mo- 
menta by replacing bq with 2sin(6g/2). Notice that in 
Eqs. ([72|, (74), and (76), the dynamic coupling modi- 



affects the excitation frequency of the magnetic superlat- 
tice. Crystal and shape anisotropies on top of the simple 
effective fields assumed above might become important 
in real structures, and their inclusion is straightforw ard . 

Let us now compare the damping {bq)^a' in E q. (74) 
with a{q) = {a±/S)q'^ corresponding to Eq. (33), which 
is the analogous quantity for the bulk. Keeping only 
the mixing conductance contribution to Eq. ( 73 1 and 
approximatin^ « p^p/in in terms of the character- 
istic Fermi momentum in the normal métal, we have 
for the g-dependent part of the damping: 



a{q) = [bqfa' 



Sd 



(77) 



up to a numerical constant. At the same time, the bulk 



a{q), corresponding to Eq. (33), can be written as 



a{q) 



(Asc/Af)^ 



(78) 



which establishes a loose formai correspondence between 
the two results. Here, l — vpT is the mean free path, A^ 
the Fermi wavelength, and the ferromagnetic cohérence 
length Asc was defined in Eq. (71 ). 



Comparing Eqs. (77l and (78), we interpret the length 
scale b <-*■ Age to describe the longest distance over which 
ferromagnetic régions can communicate via spin trans- 
fer. The length scale d ^ X characterizes momentum 
scattering relevant for spin transfer, which in the case 
of the superlattice with sharp interfaces corresponds to 
the magnetic film width d: Approximating « p|,/27r 
above, we effective ly took the normal spacers to be bal- 
listic and, because of Eq. (71), the spin transfer does 



not penetrate deep into the ferromagnetic layers, mak- 
ing possible disorder scattering there irrelevant for our 
problem. 



VI. DISCUSSION AND OUTLOOK 

Estimating the numerical value of the dimensionless q^ 
damping, according to Eq. (28l, 



a{q) = 



(Tj_q 
S 



pp/q J l + (rj_Axc)2 



(79) 



we can see that it will most likely be at most compa- 
rable or smaller than the typical g = Gilbert da mp- 



ing a ^ 10 ^, in metallic ferromagnets. Damping (79) 



fies the damping similarly to the way the static coupling 



may, however, become dominant in weak ferromagnets, 
such as diluted magnetic semiconductors. We are not 
aware of systematic expérimental investigations of the 
q^ damping in metallic ferromagnets. q^ scaling of rel- 
ative linewidth was reported in Réf. [501 for the iron-rich 
amorphous Fe9o-2;Nia;Zrio alloys. However, we are not 
certain whether the strong damping observed there can 
be attributed to the mechanism discussed in our paper. 

Another intriguing context where the physics discussed 
here can play out to be important is the current-driven 
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nonlinear ferromagnetic dynamics in mesoscopic as well 
as bulk magnetic Systems. The magnetic damping 



described by Eq. (15l can be physically thought of thc 



viscous-like spin transfer between magnetic régions pre- 
cessing slightly out-of-phase. The obvions conséquence 
of this is the enhanced damping of the inhomogeneous 
dynamics and thus the synchronization of collective mag- 
netic precession. This phenomenon was predicted in 
Réf. [7] and unambiguously observed in Ref.HTJ in the case 
of the coupled dynamics of a magnetic bilayer: When the 
two layers are tuned to similar résonance conditions, only 
the symmetric mode corresponding to the synchronized 
dynamics produces a strong response, while the antisym- 
metric mode is strongly suppressed. It is thus natural to 
suggest that the viscous magnetic damping in the con- 
tinuum limit may have far-reaching conséquences for the 
current-driven nonlinear power spectrum as that mea- 
sured in Réf. [221 This needs a further investigation. 

The rôle of electron-electron interactions was mani- 
fested in our theory through the spin Coulomb drag, 
which enhances the effective transverse spin scattering 
rate (29l. This becomes particularly important, in com- 



parison to the disorder contribution to the transverse spin 
scattering, in the limit of weak magnetic polarization.8 



We finally emphasize that the study in this paper was 
limited exclusively to weak linearized perturbations of 
the magnetic order with respect to a uniform equilib- 
rium state. When the equilibrium or out-of-equilibrium 
magnetic state is macroscopically nonuniform, as is the 
case with, e.g., the magnetic spin spirals, domain walls, 
vortices, and other topological states, the longitudinal as 
well as transverse spin currents become relevant for the 
magnetic dynamics. The longitudinal spin currents lead 
to additional contributions to the spin-transfer torques, 
modifying the magnetic équation of motion. Such spin 
torques leading to the dissipative q^ damping terms were 
discussed in Réf. 13, Thèse latter contributions to the 
magnetic damping are likely to dominate in strongly- 
textured magnetic Systems. 
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It is natural to also wonder about a possible additional 
spin current of the form ji oc didtm, which does not 
break time-reversal symmetry. Such spin current leads to a 
wave-vector-dependent correction to the effective gyromag- 
netic ratio, which is very small in practice. It parallels the 
structure of the spin pumping in magnetic nanostructures, 
which consists of the dominant dissipative pièce of the form 
m X dtm and a smaller pièce of the form dtm. The latter 
merely causes a slight rescaling of the gyromagnetic ratio. 
While the dissipative pièce of the spin pumping has been 
unambiguously established in a number of experiments,!^ 
the small correction to the gyromagnetic ratio is yet to be 
observed. 

We need to remark here that the above gauge transfor- 
mation does not affect the transverse spin current in the 
linearized theory. 



